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Abstract 

By means of self-consistent semi-empirical LCAO calculations we 
study the itinerant magnetism of (Feo.65Nio.35)i-j,Mnj, alloys for y 
between and 0.22 at T = K, neglecting only the transverse spin 
components. We find that the magnetic behaviour is quite complicated 
on a local scale. In addition to ferromagnetic behaviour, also metastable 
spin-glass-like configurations are found. In the same approach, using 
a direct numerical calculation by the Kubo- Formalism without any fit 
parameters, we also calculate the electrical conductance in the magnetic 
state and find that the y-dependence observed in the experiments is 
well reproduced by our calculations, except of an overall factor of « 5, 
by which our resistivities are too large. 



1 Introduction 



In two recent papers, ^ , we have reported progress in the calculation 
of the itinerant magnetism of ordered or disordered metallic alloys 
in a self-consistent semi-empirical LCAO approach ( LCAO= "Linear 
Combination of Atomic Orbitals"). In the first of those papers, fll], fee 
iron and particularly disordered fee Fei_a;Nia; alloys have been treated, 
with particular emphasis on the so-called Invar region around x = 0.35, 

*based on the Diploma thesis of A. Paintner, Regensburg 1994 
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where the local Fe moments respond quite sensitively on changes of 
the local neighbourhood and on the local volume. In fact, the Invar 
behaviour is essentially a consequence of strong magneto- volume effects 
(see Wassermann, for a review), which happens both in face-centered 
and also in amorphous Fe, see [|I| and |Q, respectively. 

In the second paper, we concentrated on elemental Mn and 
various Mn alloys with Ni or Fe, respectively. Here it turned out that the 
magnetism of Mn is very complicated, e.g. a-Mn has a primitive cubic 
elementary cell with 58 atoms and a very complicated antiferromagnetic 
spin configuration. This spin configuration is not even exactly known, 
but seems to be non-coUinear with at least eight non-equivalent magnetic 
sites in the elementary cell, where not only the Mn moment directions, 
but also the magnitudes of the moments are strongly different for these 
sites (see p[ for details). However, in spite of the difficulties, we finally 
obtained in a consistent description of the parameters involved in 
our semi-empirical model. Here we used 9 orbitals per atom, namely 
five 3d-, three 4p- and one 4s-orbitals, and parametrized the hopping 
elements appearing in the Hamiltonian according to a Slater-Koster 
approach with Two-Center-Integrals up to third-nearest neighbours as 
given in the book of Papaconstantopoulos, Furthermore, in our 
semi-empirical LCAO approach the atomic charges and moments are 
obtained selfconsistently from a single intra-d-orbital Hubbard integral 
U, see 1^, ^. That is, the mechanism leading to the itinerant magnetism 
is the interplay of Coulomb interaction and Pauli principle, according 
to which only two electrons with different spins can occupy the same 
orbital, but if they do, this costs an additional Coulomb repulsion U, 
so that energy can be gained by having different occupation numbers 
in the up- and down spin states, respectively. 

At the end of paper p|, we finally also obtained consistency of the 
band shift parameters for the different kinds of atoms in binary alloys 
of Mn, Ni and Fe, respectively. 

Thus we are now able to treat also the disordered ternary alloy 
(Fei_a;Nia,)i_j^Mny. Here we concentrate on the case x = 0.35, where 
without Mn one is in the above-mentioned Invar region, and where 
according to the experimentalists, with increasing Mn concentration y 
a transition to spin glass behaviour should happen (see 0, 0, ||, 3)- 
The problem of this transition is of special concern to us in this paper. 



2 



With our LCAO-formalism, we have also been able to calculate 
directly the resistivity of the alloys in the magnetic state, which has 
been measured by the same experimental group, |]IU|. 

Our paper is organized as follows: In chapter 2, we describe our 
formalism and details of the computation. In chapter 3, respectively, the 
results for the magnetic behaviour and for the resistivity are presented, 
and finally in chapter 4, we give our conclusion and discuss remaining 
open questions. 

2 Formalism 

Our alloys are modelled by 108 atoms, arranged on the sites of an 
fee cube with periodic boundary conditions. The occupation of the 
sites with Fe, Ni and Mn atoms is performed randomly according 
to the formula (Fei_a;Nia;)i_yMny. For every pair of concentrations x 
and y considered in the results below, typically ten different random 
samples have been generated, from which averages have been gained. 
Concerning the lattice constants, we took the results presented on p. 91 
of 1^, according to which the ^/-dependence of the lattice constant a can 
roughly be described by a parabola passing through a=3.59 A at y = 
and y = 0.19, with the minimum value of a = 3.584 A at y = 0.095 
in-between. 

The quantum-mechanical self-consistency equations, which have to 
be solved in our LCAO approximation, are 

m,l3 

Here / and m count the 108 sites, a and j3 the 9 orbitals per 
site, and a = ±1 the two possible orientations of the spin; for every 
spin direction, the index u counts the 9x108 possible orthonormal 
single-particle eigenvectors cj'^^ with real components c^'^^^ and the 
corresponding single-particle energies e^^\ The parameter Ui^a is the 
above-mentioned Hubbard-type intra-orbital Coulomb integral (see below), 
while ni^a,cr and n^^'^'^'* are the expectation values of the occupation 
numbers in the magnetic and paramagnetic state, respectively, and 
determined by the self-consistency equations 
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Uf{a) 

ni,.,.= E (2) 
i/=i 

In eqn. (|), the sum u = 1, i//((t) is over the occupied single- 
particle energies of the given spin-direction. In the present disordered 
magnetic system , these states are non-degenerate. 

Of course it should be mentioned that this treatment of the Hubbard 
interaction corresponds to a mean-field approximation, so that in the 
calculation of the total energy one would have to subtract the well- 
known double counting correction terms from the sum of the single 
particle energies (see e.g. @, IH); however, as already mentioned above, 
the energies calculated from our approach are in any case not very 
reliable. 

Finally, there appear parameters A/ in eqn. (|l]), which play the role 
of effective fields serving for the initialization of the magnetic state, 
if one starts from a paramagnetic state. These initializing fields are 
switched off after the first few iterations; their position dependence 
becomes an important point, when one wants to generate spin glass 
configurations (see below). 

Concerning the values of Ui^a, we assume as usual that they are 
different from only for d-states and actually depend only on the kind 
of atom considered, namely Ui^a= 5.3 eV for Fe and Ni, see |]T|, and 
2.9 eV for Mn, see §. 

The if"^-matrix in eqn. (0) describes the LCAO Hamiltonian for 
the paramagnetic system. Here, the elements with I = m, but a (3, 
vanish, while for the terms with a = f3,i. e. the self-energies, we use the 
numbers given in the book of Papaconstantopoulos, p. 289 ff. And 
as already mentioned in @], although the Fe values are taken from bcc 
Fe, they work even for fee alloys of Fe, Ni and/or Mn, if one introduces 
additional self-energy -shifts dpe^Mn = +0.9 eV and dpe^m = +0.95 
eV for the if;""-levels of Mn resp. Ni relative to the Fe levels. Of course, 
finally in the magnetic state the ni^a,<7 are in any case self-consistently 
determined through the ?7-terms in eqn. (|l]). 

Concerning the hopping matrix elements H^^, i.e. with I ^ m, we 
have again relied on the parametrization with two-center integrals and 
orthonormal orbitals in f^, and again, as in we take into account 

the fact that the distances between the sites in our alloys are different 
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from those in by modifying the two-center integrals /«,/?; 7 (e.g. with 
7 = cr, TT or 5) according to the ansatz 

Here /c = 5 for (a, P)= (d,d); k = 7/2 for (a, /5)= (p,d) or (s,d), and k = 
2 for (q!,/?)= (s,s), (s,p) or (p,p); see Harrison |jTT|. The //-parameters 
are fitted from the /-values given in for the distances there; they 
are different for the first, second and third neighbour shell. For Fe, the 
situation is somewhat more complicated: There, the 77-values for the 
three neighbour shells in fee Fe are obtained from those contained in 
the results for bcc Fe in with the interpolation procedure given in 
the appendix of . 



This determines finally the paramagnetic part of our Hamiltonian: 
For the magnetic properties, in which we are finally interested, and also 
for the charge distribution in the magnetic state, self-consistency is then 
demanded, expressed by the terms proportional to Ui^a in eqn. (|I]). It 
is important that this self-consistency does not refer only to global 
properties, but actually it refers to the local values of all charge- and 
spin densities. 

At the end, we comment on the parameter A; in eqn. ([I|): If we 
want to generate a globally ferromagnetic state, we start everywhere 
with positive values A/ = 1 eV, which are then switched off after the 
first iteration. In contrast, for the generation of spin glass states we 
use a more complicated procedure: At first, we start with a A;, which 
is 7^ only for one Fe site (for which we usually take one surrounded 
by a large number of Ni atoms, since there seems to be a rather strong 
and positive exchange interaction between Fe and Ni). Then, in five 
iterations of eqn. (|l]), a nontrivial spin configuration develops in the 
whole sample, which is already almost self-consistent concerning the 
final moment directions, but not concerning the amplitudes; i.e. the 
spin configuration obtained has different signs of the spin polarization 
at different atoms. This distribution of signs may look random: Actually, 
it is not, since it reflects correlations due to the interactions in the 
system. So we now continue in iteration number 6 with transient A; 
values of magnitude 1 eV everywhere, but with different signs at different 
/, namely the signs of the local spin polarization determined in the first 
five iterations. For the following iterations, the A/ are then switched off 
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completely, and the equations are iterated to self-consistency, which can 
take 100 iterations or even more. In this way, our spin glass configurations 
are generated, which still have self-consistently determined magnetic 
short-range order. 

3 Results 

3.1 Magnetic configuration 

In Fig. 1, we present the magnetic phase diagram as determined ex- 
perimentally by p, 1^. Here, for 0.04 < y < 0.14, one observes re- 
entrant spin glass behaviour, i.e. a transition from paramagnetic via 
ferromagnetic to spin glass behaviour, if the temperature is decreased. 
However, already at this point one should note that there exists a 
particular region between the ferromagnetic and the spin-glass regions, 
where both states apparently co-exist as metastable states. In fact, in 
Fig. 2 we show the results of our numerical simulation with ferromagnetic 
preparation: Here, up to y ~ 0.07, a nice agreement between the 
experimentally determined average moment and our results is found, 
whereas for larger y our calculated behaviour apparently continues as 
before in the ferromagnetic state with gradually decreasing average, 
whereas experimentally, there is a more drastic decrease to the spin 
glass phase. However, since in our calculation we use periodic boundary 
conditions with a periodicity length of only three fee lattice constants 
and in view of the other approximations of our method, this discrepancy 
is perhaps not astonishing. Namely, throughout the same concentration 
region, by the above-mentioned spin glass preparation, we can also 
get the already mentioned metastable spin glass states with essentially 
vanishing average moment. A typical plot of the moment distribution in 
such a spin glass state is presented in Fig. 3 for a disordered Feo.65Nio.35 
system. 

Comparing the energy of the ferromagnetic and the spin glass states, 
we find that our ferromagnetic states are somewhat more favourable 
by energies of ~ 6 meV per atom. However, as mentioned above, the 
accuracy of our method is not reliable enough to make this a significant 
statement. In fact, experimentally, see Fig. 1, the spin glass state seems 
to have lower energy; but at least the energy difference, corresponding 
to T ^ 60K, is of the same order. 
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In Fig. 4, we present histograms for the magnetic moment distributions 
of Fe-, Ni-, and Mn-atoms for 11 different Mn concentrations y of the 
ferromagnetic states calculated by our method for 
(Feo.65Nio.35)i-j/Mny. Fig. 5 exhibits the average moments versus concentration. 
Interestingly, the Ni moments have apparently always well defined 
moments around 0.5 /i^, whereas the Fe moment distributions are 
rather broad, e.g. they are spread between ~ 0.5 and ~ 2.5 /i^ , for 
y = 0.04. Of course, details of the Fe distribution depend on the Mn 
concentration y, e.g. for y = only Fe moments above 1.5 hb are 
observed, whereas for y ^ 0.05 also moments around 0.5 fiB appear, 
but generally, the ^/-dependence of the Fe-moment distributions does 
not look drastic. 

With the Mn moments itself, this is different: Some of the Mn atoms 
are polarized parallel to the majority spin direction, but the majority 
of the Mn atoms seems to be polarized antiparallel, as one can see 
from Fig. 5. Moreover, the magnitude of the Mn moments can vary 
drastically from site to site, which is not astonishing, since this was 
already true for a-Mn, [Q. 

For the spin glass configurations it is characteristic, see Fig. 6, that 
the Ni moments are no longer spread around 0.5 fiB, but around zero. 
Further, for y = 0.2 both the Fe and the Mn moments are almost 
homogeneously distributed between ^ —2 fiB and ~ +2 /i^. (For y =0, 
the Fe moments are distributed between —2.5 and +2.5 /is-) 

In Fig. 7-9 we compare our calculated distributions for the Fe mo- 
ments with hyperfine field measurements obtained by Mossbauer spec- 
troscopy experiments |13|. For y=0 and y=0.039 (Fig. 7 and 8) there is 
good agreement between the experimental hyperfine field distribution 
and the magnetic moment distribution in the ferromagnetic state. For 
higher Mn-concentrations, e.g. y=0.102, the spin glass state corresponds 
better to the experimental result (Fig. 9). For each of this three con- 
centrations we assume the same conversion factor of 14.5 T per /i^. 

According to ||13| the oscillatory structure of the hyperfine field dis- 
tribution is of no physical relevance, but is simply a consequence of the 
fit procedure used by the authors. 

A natural question is, whether the magnetic moments depend on the 
number of neighbouring atoms of the same or different kind. Therefore, 
in Fig. 10, we present some kind of plot, which allows answers on this 



7 



question: E.g., from this figure it is obvious that high Fe moments are 
favoured by a high number of Ni neighbours, and disfavoured by a 
high number of Mn neighbours, whereas the Ni moments themselves 
are rather insensitive on the neighbourhood. However, the spread of 
the distributions in Fig. 10 around the average data points should be 
kept in mind, when stating these tendencies, which can of course be 
quantified only on average, see the open circles in Fig. 10. Finally, 
it should be noted again that the magnitude of the Mn moments 
range between and ~ 2 yU^, with typical values around I fis- Fig. 11 
presents similar plots for the spin glass configurations, where the width 
of the distributions is impressive, and the tendencies concerning Ni 
neighbours in case of Fe atoms are similar. 

We have tried similar plots as in Fig. 10 and Fig. 11 also for the 2nd- 
and 3rd-nearest neighbour shell, however with non-conclusive results. 
Therefore, those plots are omitted. 

3.2 Simplified spin model 

In the following we try to describe the results for the magnetic behaviour 
with a simple molecular-field type ansatz. For the determination of the 
magnetic moments /if at the different sites / of the sample we have 
tried the equation 

42 
k=l 

In eqn. (^, the sum is over the 42 neighbours of site / up to the third- 
nearest neighbour shell, and there are nine coefficients per shell, namely 

CA{l),A{k) G {cFe,Fe', CFe,Ni', CFe,Mn', CMn,Fe', CMn,Nu CMn,Mn}y depending 

on the occupation of the sites / and k by the respective kind of atoms. 
Qualitatively, the terms CA{i),Aik) ■ /^fc can be understood as composed 
of a local susceptibility xi times a local effective field hi^k oc fi^. In 
fact, fitting the coefficients to our results for /if, which we obtained 
by our itinerant approach, with the ansatz of eqn. (^), we obtain the 
results presented in Fig. 12. Although these results depend on the 
Mn concentration y, and although the statistical error bars, which 
we have omitted to show, are rather large, namely typically ±0.06, 
some conclusions can be drawn from Fig. 12: In the ferromagnetic 
state, the magnetization of Fe atoms is strongly influenced by the Ni 
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neighbors in the first neighbour shell, and up to ?/ ~ 0.1 also by those 
of the second neighbour shell. In fact, while CFe,Ni is as large as 0.3 
for y=0 and increases up to 0.4 for larger y, Cpe^Fe is only around 

0. 06 for ?/=0 and even as low as 0.025 for y > 0.1; however, it should 
of course be taken into account that /iTVi ~ 0.5, while can be 
three or five times as large, so Ni is dominating, but not strongly. 
Moreover, the coefficient cpcMn is negative over a considerable y-range 
and particularly strong in the 2nd-neighbour shell for small y-values, 
similar as CMn,Mn, whereas CFe,Fe is positive in the ferromagnetic state, 
as mentioned above, contrary to a formerly wide-spread assumption 
([Q). Also interesting is the strong scatter of the coefficients for Mn, 

1. e. the third column. This might show that the Mn moments tend to 
get "frustrated", however a definite statement cannot be made in view 
of the small number of Mn atoms. 

Most important, however, is the fact that for the spin glass states 
(for which we omit the corresponding plots |T^) the coefficients are 



different from those of the ferromagnetic case presented in Fig. 12: The 
main difference concerns cpe^Ni-, which becomes as high as 2.3 for ?/=0, 
1.5 for ?/=0.1 and 1.4 for y = 0.2. This is understandable, since in the 
spin glass state also fim decreases in magnitude from 0.5 /ig to much 
smaller values (see above); however due to the increase of CFe,Ni, the Ni 
moments remain important for the Fe moments. A second important 
difference between the ferromagnetic and the spin glass states is that 
CFe,Fei which was positive for ferromagnetic preparation, is negative 
(^ —0.1) for spin glass preparation. So the implicit assumption made 
with the above-mentioned ansatz, namely that the CA(i),A(k) coefficients 
should not depend on the global magnetic state, is obviously not true. 
This means, one should better stay with the itinerant description, 
realizing that here the spin configuration is not fixed by given interaction 
coefficients as in Heisenberg spin glasses: Instead, the effective interactions 
depend on the global itinerant state and are different for a ferromagnetic 
and an itinerant spin-glass configuration of the system. Thus the itinerant 
description is really necessary. 

3.3 Electrical resistivity in the magnetic state 

It is perhaps not widely known that in a disordered crystalline or 
amorphous metallic system not only single particle properties, as just 
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discussed, but even two-particle properties as the electrical resistivity 
can be calculated by our semi-empirical LCAO formalism, namely by 
means of the Kubo formalism, In fact, in [|I^], the Kubo formalism 
has been successfully applied for the numerical calculation of the resistivity 
in the magnetic state of amorphous Fe/Zr alloys, and there we have 
also described how the conductivity a^^x at T=0 K as a function of the 
Fermi-energy Ej can be calculated from the formula 



Here the sum is over all u with e^'^^ G AE; further, is the volume 
per atom, c^'^^ an eigenstate of eqn. ([^) with energy e^*^-* in the interval 
AE around Ef, and rj is an effective line broadening parameter, which 
should be of the order of 1/2 the typical energy distance between 
neighbouring eigenvalues. Of course, the results depend on the choice 
of AE and of t], and additionally one has to perform sample averages 
and averages over the directions x, y, and z of the conductivity. Then, 
as we see in the Fig. 13 below, and as we have already found in |TB 
the result is essentially independent on AE and r] : 



In Fig. 13, we compare our results with the experiments in |]T0|] and 
find good agreement, concerning the dependence on the Mn concentration 
y, to the accuracy considered, over the whole parameter range. However, 
there is one discrepancy, which we do not yet understand at present, 
namely our resistivities are always too large by a factor ps; 5, compared 
with the experiment. But in view of the approximations involved in 
the evaluation of eq. @ and in view of the large scale of resistivities 
of different metals this discrepancy should perhaps not be taken too 
serious. In fact with rj ^ 0.02 eV we induce effectively inelastic scattering 
events corresponding to T ^ 200 K, and the effect of AE ^ 0.2 eV may 
even be larger, whereas the experimental values in our comparison are 



at 4 K. However, at 300 K, according to [|T^, the experimental values 



would only be enhanced with respect to those at 4 K by factors between 
1.8 (for y=0) and 1.1 for y = 0.084. Additionally, the conductivity 
should be influenced by transverse components of magnetic moments, 
which we have neglected. 
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4 Discussion 



We have studied the magnetic behaviour of Fe-Ni-Mn alloys 
in the sensitive regime around the Fe-Ni Invar- region, i.e. 
(Feo.65Nio.35)i-j/ Muj^, where with increasing Mn concentration y a transition 
from ferromagnetic to spin glass behaviour has been observed experimentally, 

Actually, with our semi-empirical LCAO approach we have found 
that both ferromagnetic states and also itinerant spin glass states can 
be prepared and are metastable in the whole concentration region 
considered. These different states have roughly the same energy, i.e. according 
to our calculation, the ferromagnetic state appears to be slightly favoured 
by ~ 6 meV per atom; however our method is not accurate enough 
for a definite statement. In any case this implies that the itinerant 
magnetism of these alloys is quite complicated and cannot simply be 
described by the Heisenberg model, since one has e.g. a broad distribution 
of different values of the magnitudes of the local atomic moments, both 
for Fe and Mn. 

Starting our iterations with a homogeneous ferromagnetic state, we 
obtained good agreement with experimental results in the ferromagnetic 
state up to a Mn concentration of ?/ ~ 0.07. For higher concentrations, 
more realistic results (e.g. the already mentioned spin-glass states with 
vanishing total magnetic moment, but non-trivial local moments) have 
been obtained systematically from a self-consistently prepared inhomogeneous 
state with correct short range magnetic order. Thus there seems to be 
some kind of preparation dependence of the equilibrium state. After 
all, this is typical for a situation with possible spin-glass behaviour and 
not in contradiction to the experimental situation, see 0. 

Finally, we have also calculated the electrical resistivity in the magnetic 
state in our LCAO formalism by an implementation of the Kubo formalism 
(see [|16|), and we obtained good agreement with experiments apart 
from the fact that our conductivities are generally too low by a factor of 
~ 5, compared with experiments. The problem, whether the conductivity 
in the magnetic state will change with the introduction of transverse 
spin components, remains to be studied in the future. 
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Figure Captions 

Fig. 1: Magnetic phase diagram for the system (Feo es^^o 35)1-1/^^1/, from 
[6,8]. 

Fig. 2: Magnetization of the alloy system (Feo.65-^?o.35)i-2/-^% versus Mn- 
conccntration y. 

Fig. 3: Typical moment distribution of a spin glass state for a disordered 
Feo.65Nio.35 alloy. The Ni moments are scaled by a factor 10. The symbols 
© above and below an atomic position show the local spin direction 
of the atoms at the intermediary state in the 6-th iteration (sec text). 

Fig. 4: Calculated histograms for the magnetic moment distribution of Fe-, 
Ni- and Mn-atoms for Mn concentrations y in the range < |/ < 0.231 
(ferromagnetic preparation) . 

Fig. 5: Average magnetic moments for the three different alloy components 
Fe, Ni and Mn. 

Fig. 6: Calculated histograms for the magnetic moment distribution of Fe-, Ni- 
and Mn-atoms for Mn concentration y e {0; 0.102; 0.204} (spin glass 
preparation). 

Fig. 7: Calculated magnetic-moment probability histogram and measured hyperfine 
field probability distribution (dashed- line fit) of the system Feo.65A^^o.35- 
The simulation was done with ferromagnetic preparation. A conversion 
factor of i7 = 14.5 T per is assumed. 

Fig. 8: Calculated magnetic-moment probability histogram (ferromagnetic 

preparation) of the system (Feo.65A^^o.35)o.963-^^o.o37 and measured hyperfine 

field probability distribution (dashed-line fit) of the system (Feo.65iVio.35)o.96i^^^o.o39- 

Fig. 9: Calculated magnetic-moment probability histogram (spin glass preparation) 
of the system (Feo.65-/V^o.35)o.898^^o.io2 and measured hyperfine field 
probability distribution (dashed-line fit) of the system 

(Feo.65iVio.35)o.919Mno.o84- 

Fig. 10: Influence of the local neighbourhood on the magnitude of the magnetic 
moments for final states with ferromagnetic preparation. In the left 
three columns the magnitude of the Fe moments is plotted versus the 
next neighbour number of Fe-, Ni- and Mn- neighbours for six different 
Mn concentrations. In the middle three columns the same is plotted for 
Ni, and in the right three columns for Mn. The dashes represent the 
single values and the circles are the averaged values. 

Fig. 11: Influence of the local neighbourhood on the magnitude of the magnetic 
moments for final states with spin glass preparation. 
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12: The coefficients CA{i),A(k) from equation (4) arc plotted against Mn 
concentration y for states with ferromagnetic preparation. The coefficients 
of the first row are related to next neighbours. Those of the second and 
third row are related to second resp. third neighbour shells. 

13: Measured and calculated conductivities against Mn concentration for 
several combinations of the parameters rj and AE. 
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